Direct neuronal reprogramming of NDUFS4 patient cells identifies the unfolded protein response as a novel general reprogramming hurdle.

Mitochondria account for essential cellular pathways, from ATP production to nucleotide metabolism, and their deficits lead to neurological disorders and contribute to the onset of age-related diseases. Direct neuronal reprogramming aims at replacing neurons lost in such conditions, but very little is known about the impact of mitochondrial dysfunction on the direct reprogramming of human cells. Here, we explore the effects of mitochondrial dysfunction on the neuronal reprogramming of induced pluripotent stem cell (iPSC)-derived astrocytes carrying mutations in the NDUFS4 gene, important for Complex I and associated with Leigh syndrome. This led to the identification of the unfolded protein response as a major hurdle in the direct neuronal conversion of not only astrocytes and fibroblasts from patients but also control human astrocytes and fibroblasts. Its transient inhibition potently improves reprogramming by influencing the mitochondria-endoplasmic-reticulum-stress-mediated pathways. Taken together, disease modeling using patient cells unraveled novel general hurdles and ways to overcome these in human astrocyte-to-neuron reprogramming.

Leigh syndrome is the main clinical characteristic of PTCD3 deficiency.

Mitochondrial translation defects are a continuously growing group of disorders showing a large variety of clinical symptoms including a wide range of neurological abnormalities. To date, mutations in PTCD3, encoding a component of the mitochondrial ribosome, have only been reported in a single individual with clinical evidence of Leigh syndrome. Here, we describe three additional PTCD3 individuals from two unrelated families, broadening the genetic and phenotypic spectrum of this disorder, and provide definitive evidence that PTCD3 deficiency is associated with Leigh syndrome. The patients presented in the first months of life with psychomotor delay, respiratory insufficiency and feeding difficulties. The neurologic phenotype included dystonia, optic atrophy, nystagmus and tonic-clonic seizures. Brain MRI showed optic nerve atrophy and thalamic changes, consistent with Leigh syndrome. WES and RNA-seq identified compound heterozygous variants in PTCD3 in both families: c.[1453-1G>C];[1918C>G] and c.[710del];[902C>T]. The functional consequences of the identified variants were determined by a comprehensive characterization of the mitochondrial function. PTCD3 protein levels were significantly reduced in patient fibroblasts and, consistent with a mitochondrial translation defect, a severe reduction in the steady state levels of complexes I and IV subunits was detected. Accordingly, the activity of these complexes was also low, and high-resolution respirometry showed a significant decrease in the mitochondrial respiratory capacity. Functional complementation studies demonstrated the pathogenic effect of the identified variants since the expression of wild-type PTCD3 in immortalized fibroblasts restored the steady-state levels of complexes I and IV subunits as well as the mitochondrial respiratory capacity. Additionally, minigene assays demonstrated that three of the identified variants were pathogenic by altering PTCD3 mRNA processing. The fourth variant was a frameshift leading to a truncated protein. In summary, we provide evidence of PTCD3 involvement in human disease confirming that PTCD3 deficiency is definitively associated with Leigh syndrome.

Diagnostic Odyssey in an Adult Patient with Ophthalmologic Abnormalities and Hearing Loss: Contribution of RNA-Seq to the Diagnosis of a PEX1 Deficiency.

Peroxisomal biogenesis disorders (PBDs) are a heterogeneous group of genetic diseases. Multiple peroxisomal pathways are impaired, and very long chain fatty acids (VLCFA) are the first line biomarkers for the diagnosis. The clinical presentation of PBDs may range from severe, lethal multisystemic disorders to milder, late-onset disease. The vast majority of PBDs belong to Zellweger Spectrum Disordes (ZSDs) and represents a continuum of overlapping clinical symptoms, with Zellweger syndrome being the most severe and Heimler syndrome the less severe disease. Mild clinical conditions frequently present normal or slight biochemical alterations, making the diagnosis of these patients challenging. In the present study we used a combined WES and RNA-seq strategy to diagnose a patient presenting with retinal dystrophy as the main clinical symptom. Results showed the patient was compound heterozygous for mutations in PEX1. VLCFA were normal, but retrospective analysis of lysosphosphatidylcholines (LPC) containing C22:0-C26:0 species was altered. This simple test could avoid the diagnostic odyssey of patients with mild phenotype, such as the individual described here, who was diagnosed very late in adult life. We provide functional data in cell line models that may explain the mild phenotype of the patient by demonstrating the hypomorphic nature of a deep intronic variant altering PEX1 mRNA processing.

Bi-allelic LETM1 variants perturb mitochondrial ion homeostasis leading to a clinical spectrum with predominant nervous system involvement.

Leucine zipper-EF-hand containing transmembrane protein 1 (LETM1) encodes an inner mitochondrial membrane protein with an osmoregulatory function controlling mitochondrial volume and ion homeostasis. The putative association of LETM1 with a human disease was initially suggested in Wolf-Hirschhorn syndrome, a disorder that results from de novo monoallelic deletion of chromosome 4p16.3, a region encompassing LETM1. Utilizing exome sequencing and international gene-matching efforts, we have identified 18 affected individuals from 11 unrelated families harboring ultra-rare bi-allelic missense and loss-of-function LETM1 variants and clinical presentations highly suggestive of mitochondrial disease. These manifested as a spectrum of predominantly infantile-onset (14/18, 78%) and variably progressive neurological, metabolic, and dysmorphic symptoms, plus multiple organ dysfunction associated with neurodegeneration. The common features included respiratory chain complex deficiencies (100%), global developmental delay (94%), optic atrophy (83%), sensorineural hearing loss (78%), and cerebellar ataxia (78%) followed by epilepsy (67%), spasticity (53%), and myopathy (50%). Other features included bilateral cataracts (42%), cardiomyopathy (36%), and diabetes (27%). To better understand the pathogenic mechanism of the identified LETM1 variants, we performed biochemical and morphological studies on mitochondrial K+/H+ exchange activity, proteins, and shape in proband-derived fibroblasts and muscles and in Saccharomyces cerevisiae, which is an important model organism for mitochondrial osmotic regulation. Our results demonstrate that bi-allelic LETM1 variants are associated with defective mitochondrial K+ efflux, swollen mitochondrial matrix structures, and loss of important mitochondrial oxidative phosphorylation protein components, thus highlighting the implication of perturbed mitochondrial osmoregulation caused by LETM1 variants in neurological and mitochondrial pathologies.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics.

BACKGROUND: Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics. METHODS: We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage. RESULTS: We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions. CONCLUSION: Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

Genetic basis of mitochondrial diseases.

Mitochondrial disorders are monogenic disorders characterized by a defect in oxidative phosphorylation and caused by pathogenic variants in one of over 340 different genes. The implementation of whole-exome sequencing has led to a revolution in their diagnosis, duplicated the number of associated disease genes, and significantly increased the diagnosed fraction. However, the genetic etiology of a substantial fraction of patients exhibiting mitochondrial disorders remains unknown, highlighting limitations in variant detection and interpretation, which calls for improved computational and DNA sequencing methods, as well as the addition of OMICS tools. More intriguingly, this also suggests that some pathogenic variants lie outside of the protein-coding genes and that the mechanisms beyond the Mendelian inheritance and the mtDNA are of relevance. This review covers the current status of the genetic basis of mitochondrial diseases, discusses current challenges and perspectives, and explores the contribution of factors beyond the protein-coding regions and monogenic inheritance in the expansion of the genetic spectrum of disease.

Detection of aberrant splicing events in RNA-seq data using FRASER.

Aberrant splicing is a major cause of rare diseases.  However, its prediction from genome sequence alone remains in most cases inconclusive. Recently, RNA sequencing has proven to be an effective complementary avenue to detect aberrant splicing. Here, we develop FRASER, an algorithm to detect aberrant splicing from RNA sequencing data. Unlike existing methods, FRASER captures not only alternative splicing but also intron retention events. This typically doubles the number of detected aberrant events and identified a pathogenic intron retention in MCOLN1 causing mucolipidosis. FRASER automatically controls for latent confounders, which are widespread and affect sensitivity substantially. Moreover, FRASER is based on a count distribution and multiple testing correction, thus reducing the number of calls by two orders of magnitude over commonly applied z score cutoffs, with a minor loss of sensitivity. Applying FRASER to rare disease diagnostics is demonstrated by reprioritizing a pathogenic aberrant exon truncation in TAZ from a published dataset. FRASER is easy to use and freely available.

Detection of aberrant gene expression events in RNA sequencing data.

RNA sequencing (RNA-seq) has emerged as a powerful approach to discover disease-causing gene regulatory defects in individuals affected by genetically undiagnosed rare disorders. Pioneering studies have shown that RNA-seq could increase the diagnosis rates over DNA sequencing alone by 8-36%, depending on the disease entity and tissue probed. To accelerate adoption of RNA-seq by human genetics centers, detailed analysis protocols are now needed. We present a step-by-step protocol that details how to robustly detect aberrant expression levels, aberrant splicing and mono-allelic expression in RNA-seq data using dedicated statistical methods. We describe how to generate and assess quality control plots and interpret the analysis results. The protocol is based on the detection of RNA outliers pipeline (DROP), a modular computational workflow that integrates all the analysis steps, can leverage parallel computing infrastructures and generates browsable web page reports.

Corrigendum: ncRNAs: New Players in Mitochondrial Health and Disease?

[This corrects the article DOI: 10.3389/fgene.2020.00095.].

ncRNAs: New Players in Mitochondrial Health and Disease?

The regulation of mitochondrial proteome is unique in that its components have origins in both mitochondria and nucleus. With the development of OMICS technologies, emerging evidence indicates an interaction between mitochondria and nucleus based not only on the proteins but also on the non-coding RNAs (ncRNAs). It is now accepted that large parts of the non-coding genome are transcribed into various ncRNA species. Although their characterization has been a hot topic in recent years, the function of the majority remains unknown. Recently, ncRNA species microRNA (miRNA) and long-non coding RNAs (lncRNA) have been gaining attention as direct or indirect modulators of the mitochondrial proteome homeostasis. These ncRNA can impact mitochondria indirectly by affecting transcripts encoding for mitochondrial proteins in the cytoplasm. Furthermore, reports of mitochondria-localized miRNAs, termed mitomiRs, and lncRNAs directly regulating mitochondrial gene expression suggest the import of RNA to mitochondria, but also transcription from the mitochondrial genome. Interestingly, ncRNAs have been also shown to hide small open reading frames (sORFs) encoding for small functional peptides termed micropeptides, with several examples reported with a role in mitochondria. In this review, we provide a literature overview on ncRNAs and micropeptides found to be associated with mitochondrial biology in the context of both health and disease. Although reported, small study overlap and rare replications by other groups make the presence, transport, and role of ncRNA in mitochondria an attractive, but still challenging subject. Finally, we touch the topic of their potential as prognosis markers and therapeutic targets.

SSBP1 mutations cause mtDNA depletion underlying a complex optic atrophy disorder.

Inherited optic neuropathies include complex phenotypes, mostly driven by mitochondrial dysfunction. We report an optic atrophy spectrum disorder, including retinal macular dystrophy and kidney insufficiency leading to transplantation, associated with mitochondrial DNA (mtDNA) depletion without accumulation of multiple deletions. By whole-exome sequencing, we identified mutations affecting the mitochondrial single-strand binding protein (SSBP1) in 4 families with dominant and 1 with recessive inheritance. We show that SSBP1 mutations in patient-derived fibroblasts variably affect the amount of SSBP1 protein and alter multimer formation, but not the binding to ssDNA. SSBP1 mutations impaired mtDNA, nucleoids, and 7S-DNA amounts as well as mtDNA replication, affecting replisome machinery. The variable mtDNA depletion in cells was reflected in severity of mitochondrial dysfunction, including respiratory efficiency, OXPHOS subunits, and complex amount and assembly. mtDNA depletion and cytochrome c oxidase-negative cells were found ex vivo in biopsies of affected tissues, such as kidney and skeletal muscle. Reduced efficiency of mtDNA replication was also reproduced in vitro, confirming the pathogenic mechanism. Furthermore, ssbp1 suppression in zebrafish induced signs of nephropathy and reduced optic nerve size, the latter phenotype complemented by WT mRNA but not by SSBP1 mutant transcripts. This previously unrecognized disease of mtDNA maintenance implicates SSBP1 mutations as a cause of human pathology.

Bi-Allelic UQCRFS1 Variants Are Associated with Mitochondrial Complex III Deficiency, Cardiomyopathy, and Alopecia Totalis.

Isolated complex III (CIII) deficiencies are among the least frequently diagnosed mitochondrial disorders. Clinical symptoms range from isolated myopathy to severe multi-systemic disorders with early death and disability. To date, we know of pathogenic variants in genes encoding five out of 10 subunits and five out of 13 assembly factors of CIII. Here we describe rare bi-allelic variants in the gene of a catalytic subunit of CIII, UQCRFS1, which encodes the Rieske iron-sulfur protein, in two unrelated individuals. Affected children presented with low CIII activity in fibroblasts, lactic acidosis, fetal bradycardia, hypertrophic cardiomyopathy, and alopecia totalis. Studies in proband-derived fibroblasts showed a deleterious effect of the variants on UQCRFS1 protein abundance, mitochondrial import, CIII assembly, and cellular respiration. Complementation studies via lentiviral transduction and overexpression of wild-type UQCRFS1 restored mitochondrial function and rescued the cellular phenotype, confirming UQCRFS1 variants as causative for CIII deficiency. We demonstrate that mutations in UQCRFS1 can cause mitochondrial disease, and our results thereby expand the clinical and mutational spectrum of CIII deficiencies.

Homozygous loss-of-function variants of TASP1, a gene encoding an activator of the histone methyltransferases KMT2A and KMT2D, cause a syndrome of developmental delay, happy demeanor, distinctive facial features, and congenital anomalies.

We report four unrelated children with homozygous loss-of-function variants in TASP1 and an overlapping phenotype comprising developmental delay with hypotonia and microcephaly, feeding difficulties with failure-to-thrive, recurrent respiratory infections, cardiovascular malformations, cryptorchidism, happy demeanor, and distinctive facial features. Two children had a homozygous founder deletion encompassing exons 5-11 of TASP1, the third had a homozygous missense variant, c.701 C>T (p.Thr234Met), affecting the active site of the encoded enzyme, and the fourth had a homozygous nonsense variant, c.199 C>T (p.Arg67*). TASP1 encodes taspase 1 (TASP1), which is responsible for cleaving, thus activating, the lysine methyltransferases KMT2A and KMT2D, which are essential for histone methylation and transcription regulation. The consistency of the phenotype, the critical biological function of TASP1, the deleterious nature of the TASP1 variants, and the overlapping features with Wiedemann-Steiner and Kabuki syndromes respectively caused by pathogenic variants in KMT2A and KMT2D all support that TASP1 is a disease-related gene.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues.

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Bi-allelic ADPRHL2 Mutations Cause Neurodegeneration with Developmental Delay, Ataxia, and Axonal Neuropathy.

ADP-ribosylation is a reversible posttranslational modification used to regulate protein function. ADP-ribosyltransferases transfer ADP-ribose from NAD+ to the target protein, and ADP-ribosylhydrolases, such as ADPRHL2, reverse the reaction. We used exome sequencing to identify five different bi-allelic pathogenic ADPRHL2 variants in 12 individuals from 8 families affected by a neurodegenerative disorder manifesting in childhood or adolescence with key clinical features including developmental delay or regression, seizures, ataxia, and axonal (sensori-)motor neuropathy. ADPRHL2 was virtually absent in available affected individuals' fibroblasts, and cell viability was reduced upon hydrogen peroxide exposure, although it was rescued by expression of wild-type ADPRHL2 mRNA as well as treatment with a PARP1 inhibitor. Our findings suggest impaired protein ribosylation as another pathway that, if disturbed, causes neurodegenerative diseases.

OCR-Stats: Robust estimation and statistical testing of mitochondrial respiration activities using Seahorse XF Analyzer.

The accurate quantification of cellular and mitochondrial bioenergetic activity is of great interest in medicine and biology. Mitochondrial stress tests performed with Seahorse Bioscience XF Analyzers allow the estimation of different bioenergetic measures by monitoring the oxygen consumption rates (OCR) of living cells in multi-well plates. However, studies of the statistical best practices for determining aggregated OCR measurements and comparisons have been lacking. Therefore, to understand how OCR behaves across different biological samples, wells, and plates, we performed mitochondrial stress tests in 126 96-well plates involving 203 fibroblast cell lines. We show that the noise of OCR is multiplicative, that outlier data points can concern individual measurements or all measurements of a well, and that the inter-plate variation is greater than the intra-plate variation. Based on these insights, we developed a novel statistical method, OCR-Stats, that: i) robustly estimates OCR levels modeling multiplicative noise and automatically identifying outlier data points and outlier wells; and ii) performs statistical testing between samples, taking into account the different magnitudes of the between- and within-plate variations. This led to a significant reduction of the coefficient of variation across plates of basal respiration by 45% and of maximal respiration by 29%. Moreover, using positive and negative controls, we show that our statistical test outperforms the existing methods, which suffer from an excess of either false positives (within-plate methods), or false negatives (between-plate methods). Altogether, this study provides statistical good practices to support experimentalists in designing, analyzing, testing, and reporting the results of mitochondrial stress tests using this high throughput platform.

NDUFB8 Mutations Cause Mitochondrial Complex I Deficiency in Individuals with Leigh-like Encephalomyopathy.

Respiratory chain complex I deficiency is the most frequently identified biochemical defect in childhood mitochondrial diseases. Clinical symptoms range from fatal infantile lactic acidosis to Leigh syndrome and other encephalomyopathies or cardiomyopathies. To date, disease-causing variants in genes coding for 27 complex I subunits, including 7 mitochondrial DNA genes, and in 11 genes encoding complex I assembly factors have been reported. Here, we describe rare biallelic variants in NDUFB8 encoding a complex I accessory subunit revealed by whole-exome sequencing in two individuals from two families. Both presented with a progressive course of disease with encephalo(cardio)myopathic features including muscular hypotonia, cardiac hypertrophy, respiratory failure, failure to thrive, and developmental delay. Blood lactate was elevated. Neuroimaging disclosed progressive changes in the basal ganglia and either brain stem or internal capsule. Biochemical analyses showed an isolated decrease in complex I enzymatic activity in muscle and fibroblasts. Complementation studies by expression of wild-type NDUFB8 in cells from affected individuals restored mitochondrial function, confirming NDUFB8 variants as the cause of complex I deficiency. Hereby we establish NDUFB8 as a relevant gene in childhood-onset mitochondrial disease.

Biallelic Mutations in SLC1A2; an Additional Mode of Inheritance for SLC1A2-Related Epilepsy.

Recently, heterozygous de novo mutations in SCL1A2 have been reported to underlie severe early-onset epileptic encephalopathy. In one male presenting with epileptic seizures and visual impairment, we identified a novel homozygous splicing variant in SCL1A2 (c.1421 + 1G > C) by using exome sequencing. Functional studies on cDNA level confirmed a consecutive loss of function. Our findings suggest that not only de novo mutations but also biallelic variants in SLC1A2 can cause epilepsy and that there is an additional autosomal recessive mode of inheritance. These findings also contribute to the understanding of the genetic mechanism of autosomal dominant SLC1A2-related epileptic encephalopathy as they exclude haploinsufficiency as exclusive genetic mechanism.